%% contourlet
% [im,map]=imread('baboon.bmp');
close all
% [im,map]=imread('BARB.bmp');
[im,map]=imread('standard_lena.bmp');
im = double(im) / 256;
ratio1=5;
 option = 'auto' ;
  nlevels = [2, 4];    % Decomposition level
imrec1=nlademo3( im, ratio1,option,nlevels );
ratio2=10;
imrec2=nlademo3( im, ratio2,option,nlevels );
ratio3=20;
imrec3=nlademo3( im, ratio3,option,nlevels );
%%  curvelet
% [im,map]=imread('baboon.bmp');
close all
[im,map]=imread('BARB.bmp');
% [im,map]=imread('standard_lena.bmp');
ratio1=0.05;
Y=curveletdemo(im, ratio1);
ratio2=0.1;
Y=curveletdemo(im, ratio2);
ratio3=0.2;
Y=curveletdemo(im, ratio3);
%% wavelet
close all
[X,map]=imread('baboon.bmp');
A2=waveletdemo(X,0.05);
A2=waveletdemo(X,0.1);
A2=waveletdemo(X,0.2);
[X,map]=imread('BARB.bmp');
A2=waveletdemo(X,0.05);
A2=waveletdemo(X,0.1);
A2=waveletdemo(X,0.2);
[X,map]=imread('standard_lena.bmp');
A2=waveletdemo(X,0.05);
A2=waveletdemo(X,0.1);
A2=waveletdemo(X,0.2);

